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The strength of electronic correlation effects in the spin-dependent electronic structure of fer- 
romagnetic bcc Fe(llO) has been investigated by means of spin and angle-resolved photoemission 
spectroscopy. The experimental results are compared to theoretical calculations within the three- 
body scattering approximation and within the dynamical mean-field theory, together with one-step 
model calculations of the photoemission process. This comparison indicates that the present state 
of the art many-body calculations, although improving the description of correlation effects in Fe, 
give too small mass renormalizations and scattering rates thus demanding more refined many-body 
theories including non-local fluctuations. 

PACS numbers: 75.70.Rf, 79.60.Bm, 73.20.At, 71.15.Mb, 75.50.Bb 



Since more than half a century it is clear that the band- 
structure together with exchange and correlation effects 
play an important role for the appearance of ferromag- 
netism in 3(i transition metals and their alloys [1]. A 
first step toward an understanding of the electronic struc- 
ture of these metals has been achieved by calculations 
of the single-particle band dispersion [^(k)] within the 
density functional theory (DFT) in the local spin-density 
approximation (LSDA) [2 which takes into account cor- 
relation effects only in a limited extent. It soon turned 
out that for the ferromagnetic 3(i transition metals such 
as Fe, Co, and Ni, calculations beyond DFT-based theo- 
ries have to be developed to take into account many-body 
interaction, i.e., correlation effects, which normally are 
described by the energy and momentum dependent com- 
plex self-energy function I](£^,k). Here the real part 
is related to the mass enhancement while the imaginary 
part describes the scattering rate. One of the success- 
ful schemes for correlated electron systems is the dynami- 
cal mean- field theory (DMFT). It replaces the problem of 
describing correlation effects in a periodic lattice by a cor- 
related impurity coupled to a self-consistent bath [3] . An 
alternative approach is the three-body scattering (3BS) 
approximation which takes into account the scattering of 
a hole into an Auger-like excitation in the valence band, 
formed by one hole plus an electron- hole excitation [4|. 
Such many-body calculations allowed the qualitative de- 
scription of the quenching of majority-channel quasipar- 
ticle excitations in Co [5 or the narrowing of the Ni 3(i 
band [6 . While the above mentioned many-body theo- 
ries give an improved description of the electronic struc- 



ture, the central question is, whether they also lead to a 
quantitative agreement with experiments. 

Angle-resolved photoemission spectroscopy (ARPES) 
is a powerful method to determine the spectral function 
and by comparison with the bare-particle band structure 
(usually approximated by DFT band structure calcula- 
tions) to obtain the self-energy [7 . Moreover, the spin- 
resolved version of this method is very useful to disen- 
tangle the complex electronic structure of ferromagnets, 
in particular for systems with a strong overlap between 
majority and minority bands [8 . For ferromagnetic bcc 
Fe, which is the focus of the present work, several (spin- 
resolved) ARPES studies have been presented in the lit- 
erature [i nni nn na US] . 

In this letter we present (spin-resolved) ARPES spec- 
tra of ferromagnetic bcc Fe(llO). These experimental 
results are compared to theoretical many-body calcula- 
tions (DMFT and 3BS) as mentioned above, including 
a full calculation of the one-step photoemission process. 
Our main conclusion is that a quantitative agreement, 
in particular concerning the line widths, is not observed. 
This clearly demonstrates the demand for further non- 
local many-body theories. 

Experiments have been performed at room tempera- 
ture in normal emission using a hemispherical SPECS 
Phoibos 150 electron energy analyzer and undulator radi- 
ation from the UE112-PGM1 beamline at BESSY II. The 
linear polarization of the photon beam could be switched 
from horizontal to vertical, which means that a more p- 
(A||[lTO]) or s- (A||[001]) character of the incident light 
was achieved, respectively. For spin analysis, a Rice Uni- 
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versity Mott-type spin polarimeter has been used, oper- 
ated at 26 kV [14^ . The angular resolution of the equip- 
ment was 1°, and the energy resolution was 100 meV. 
The Fe(llO) surface has been prepared on W(llO) by 
deposition of 20 ML Fe and post-annealing. The struc- 
tural quality of the films has been checked by low-energy 
electron diffraction. The Fe(llO) film was remanently 
magnetized in the film plane along the [110] easy axis. 

Theoretical calculations of the ARPES spectra have 
been performed within the framework of the Korringa- 
Kohn-Rostoker multiple scattering theory [15 . Spin- 
orbit coupling and exchange splitting were treated on 
equal footing in a fully relativistic theory. To account 
for electronic correlations beyond the LSDA approxima- 
tion, a site-diagonal, local and complex energy-dependent 
self-energy T^dmft was introduced self-consistently [16]. 
In most of the calculations presented for the comparison 
with the ARPES results, unless specified, we use for the 
averaged on-site Coulomb interaction U a value [7=1.5 
eV which is within the experimental value U ^1 eV [17j 
and a value U ~2 eV derived from theoretical studies 
[ISl [19] . It is usually accepted that the averaged on-site 
exchange interaction J coincides with its atomic value 
J -0.9 eV [20 . Using this LSDA+DMFT approach the 
spectral function A{E, k) could be calculated for a semi- 
infinite lattice structure. For a quantitative comparison 
with the measured ARPES spectra, it is inevitable to 
take into account the wave- vector and energy dependent 
transition matrix elements calculated within the one-step 
model of photoemission (ISM) [21 . This describes the 
excitation process, the transport of the photoelectrons to 
the surface as well as the escape into the vacuum. For a 
quantitative description of surface states and resonances 
a realistic approach for the surface barrier is essential. Fi- 
nally, analogous calculations of the ARPES signal have 
been performed on the basis of the 3BS self-energy func- 
tion. The two calculations (3BS and DMFT) do not show 
big differences. We emphasize that in all the calculations 
the broadening of the peaks due to final state effects is 
included. 

Figures 1(a) and (b) display a comparison be- 
tween spin-integrated ARPES data and theoretical 
LSDA+DMFT+ISM calculations of Fe(llO) along the 
FN direction of the bulk Brillouin zone (BZ) with p- 
polarization. The k values were calculated from the used 
photon energies ranging from 25eV to 100 eV, using an 
inner potential V=14.5 eV. Figure 2 displays analogous 
spin-resolved ARPES data for p and s-polarized photons 
together with LSDA+DMFT+ISM calculations. 

Near the F point (k~0.06 FN), the intense peak close to 
the Fermi level corresponds to a 3 minority surface res- 
onance. Experimentally, its S3 bulk component crosses 
the Fermi level at k ~0.33 FN, leading to a reversal of 
the measured spin-polarization and to a strong reduction 
of the intensity at k =0.68 FN in the minority channel. 
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FIG. 1: (a) Experimental spin-integrated photoemission spec- 
tra of the Fe(llO) surface measured with p-polarization in 
normal emission along the FN direction of the bulk Brillouin 
zone. The curves are labeled by the wave vectors in units 
of rN=1.55 A~^. (b) Corresponding calculations obtained by 
the LSDA+DMFT+ISM method which includes correlations, 
matrix elements and surface effects. 

in agreement with the theoretical results [Fig. 2(b) and 
(d)]. The peak at the binding energy BE~0.7 eV, visible 
mainly for p-polarization in a large range of wave vectors 
between F and N can be assigned to almost degener- 
ate 4 bulk-like majority states [Figs. 1 and 2(a) and 

(c)]. For s-polarization [Fig. 2(b) and (d)], a E3 feature 
at BE^l.l eV dominates the spectrum at the F point. 
For p-polarization its degenerate states form a shoul- 
der around the same BE. The broad feature around 2.2 
eV, visible at various k points, but not at the N point, 
is related to a majority 3 surface state (see below). 
Around the N-point (0.76< k <1.0) and at BE>3 eV 
[Figs. 1(a) and 2(a)] we observe a band having strong 
sp character. The pronounced difference between its the- 
oretical and experimental intensity distributions can be 
attributed to the fact that in the present calculations only 
local Coulomb repulsion between d electrons is consid- 
ered, without additional lifetime effects for the sp bands. 
Finally, we notice that the background intensity of the 
spectrum at k=0.66 FN, corresponding to a photon en- 
ergy of 55 eV, is strongly increasing by the appearance 
of the Fe 3p resonance. 

Comparing the experimental results from spin- 
integrated and spin-resolved ARPES measurements with 
LSDA+DMFT+ISM results, we obtain at low BE good 
agreement for many of the peak positions. This is also 
demonstrated in Fig. 3(a) and (b) were we compare 
the experimental peak positions with the LSDA+DMFT 
spectral function. Similar calculations based on the 
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FIG. 2: (color online). Analogous data as in Fig. 1(a) and 
(b) but now spin-resolved, (a), (b): Experiment [black (red) 
upwards (downwards) triangles for majority (minority) elec- 
tronic states], (c), (d): LSDA+DMFT+ISM theory [black 
(dark) and red (light) lines for majority and minority elec- 
trons, respectively]. (a),(c) for p- and (b), (d) for s- polariza- 
tion. The symbols used for peak assignment are the same as 
in Fig. 1. 



LSDA+3BS scheme are compared with the experimental 
data in Fig. 3(c) and (d). Since the theoretical calcula- 
tions do not show big differences, also the LSDA+3BS 
spectral function agrees well at low BE with the experi- 
mental peak positions. 

On the other hand, quantitative agreement cannot be 
achieved for higher BE. In particular, the calculated spec- 
tral weight near F for the 3 bands is in between the 
experimental features at 1.2 eV and 2.2 eV. Assuming 
negligible correlation effects would move the calculated 
feature to the LDA value at BE=2.2 eV. Thus the ex- 
perimental peak at 2.2 eV could be assigned to the bulk 
T\ 3 bands. However, a complete neglect of correlation 
effects in Fe would make the overall comparison between 
theory and experiment much worse. Thus we interpret 
the experimental peak at BE=2.2 eV by a 3 surface 
state in agreement with previous experimental and the- 
oretical studies [22]. Our theoretical results confirm this 
view since we clearly observe how changes in the surface 
barrier potential induce additional shifts in its BE posi- 
tion. Thus from the data shown in Fig. 3 we conclude 
that correlation effects in the present calculations using 
[7=1.5 eV are underestimated and that a stronger band 
narrowing is needed to achieve agreement between theory 
and experiment. 

Figure 4 demonstrates this for data close to the F 
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FIG. 3: (color online). Spectral functions of Fe(llO) and 
photoemission peak positions obtained from the spin-resolved 
measurements for different polarizations (O for horizon- 
tal and O for vertical polarization). Results obtained by 
LSDA+DMFT [(a),(b)] and by LSDA+3BS [(c), (d)] methods 
for majority and minority electronic states, respectively. 



point (/c=0.05 FN). Figure 4(a) compares the experimen- 
tal spin-integrated ARPES spectra with LSDA calcula- 
tions broadened with the experimental energy resolution, 
a LSDA+DMFT calculation, and a LSDA+DMFT+ISM 
calculation. At low BE, perfect agreement between the- 
ory and experiment is achieved for the minority 3 sur- 
face resonance. For the bulk Y\ 4 peak at BE=0.65 eV, 
the agreement is less satisfactory. This holds also for the 
3 peak which appears in LSDA at 2.2 eV. Due to cor- 
relation effects it is shifted in the experiment to BE?^ 
1.2 eV causing in Fig. 4(a) at that energy a shoulder 
and in Fig. 2(b) for k =0.06 FN a peak. The difference 
between the ARPES data and LSDA calculations can 
be explained with a linear 5Ri;=0.7^ corresponding to 
a mass enhancement m* /mo =1.7, where mo is the bare 
particle mass. This experimental mass renormalization 
in the energy range BE< 1.8 eV should be contrasted 
to m * /mo ~1.25 derived for [7=1.5 eV in the present 
work. One may speculate that the difference between 
the ARPES and the LSDA+DMFT+ISM peak positions 
could be reduced by choosing a higher U (e.g. U ^ A eV). 
This value, however, is outside the previously mentioned 
range and in addition we still are left with the problem 
that the calculated width is far too small compared with 
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FIG. 4: (a) Comparison between the spin-integrated 
experimental spectra at the F point for p-polarization 
with the single-particle LSDA-based calculation, the 
LSDA+DMFT spectra, and the LSDA+DMFT+ISM spec- 
tra.(b) LSDA+DMFT+ISM calculations for U=1.5 to 3 eV. 



the experimental value. It is remarkable that the width 
does not increase with increasing U [Fig. 4(b)]. This can 
be explained in terms of an energy dependent and by 
the fact that with increasing U the peak position moves 
to lower BE leading to an almost constant ^S. The en- 
ergy dependence of also leads to slightly asymmetric 
peaks in the spectral function (Fig. 4(b)). 

We have observed this additional broadening not only 
at the F point but also at other k values. We have fitted 
the spin-resolved energy distribution curves also at vari- 
ous k values by a sum of Lorentzians plus a background. 
Although we are aware that such an evaluation is prob- 
lematic due to the formation of asymmetric Lorentzians 
discussed above, we have tried to extract as a func- 
tion of BE and k (not shown). From this evaluation 
we obtain in a first approximation a /c-independent ^E, 
which is roughly a factor two bigger than the calculated 

^^DMFT- 

We point out that this additional broadening cannot 
be caused by final state effects [23 , since those are fully 
taken into account in the ISM calculations. Furthermore, 
since the 3d bands show a rather flat dispersion, in partic- 
ular close to the high symmetry points points of the BZ, 
the final state broadening can almost be neglected when 
compared with the experimental resolution. Moreover, 
the final state effects would cause a broadening which 
for a small initial state dispersion, would be constant as 
a function of the BE. In this context, we mention that 
the broadening due to defect scattering is believed to be 



energy independent, as well. 

In summary, the comparison of (spin-resolved) ARPES 
data of ferromagnetic bcc Fe(llO) to state of the art 
many-body LSDA+DMFT+ISM and LSDA+3BS+1SM 
calculations indicates that these theories lead to an im- 
portant improvement compared to LSDA. However, we 
have demonstrated that these theories, at least in the 
present implementation, underestimate the mass renor- 
malization and in particular the scattering rates in Fe, 
a typical 3d multiband transition metal. This result de- 
mands for more refined many-body calculations, possibly 
with non-local schemes, which take into account non- 
local fluctuations in multiband systems. These conclu- 
sions are also important for other correlated multiband 
systems such as Fe-pnictide high- Tc superconductors. 
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